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ABSTRACT 

We consider the effects of supernovae (SNe) on accretion and star formation in a 
massive gaseous disk in a large primeval galaxy. The gaseous disk we envisage, roughly 
1 kpc in size with ^IO^Mq of gas, could have formed as a result of galaxy mergers where 
tidal interactions removed angular momentum from gas at larger radius and thereby 
concentrated it within the central ~ Ikpc region. We find that SNe lead to accretion 
in the disk at a rate of roughly 0.1-1 Mq yr~^ and induce star formation at a rate of 
~ 10-100 M.Q per year which contributes to the formation of a bulge; a part of the stellar 
velocity dispersion is due to SNa shell speed from which stars are formed and a part due 
to the repeated action of stochastic gravitational field of SNe remnant network on stars. 
The rate of SNa in the inner kpc is shown to be self regulating, and it cycles through 
phases of low and high activity. The supernova-assisted accretion transports gas from 
about one kpc to within a few pc of the center. If this accretion were to continue down 
to the central black hole then the resulting ratio of BH mass to the stellar mass in the 
bulge would be of order ~ 10^^-10"'^, in line with the observed Magorrian relation. 

Subject headings: accretion: theory, method: analytical - supernovae - stars: formation 
- galaxies: bulges 

1. Introduction 

The CO observations of ultra-luminous infra-red galaxies (ULIGals) find the gas mass in the 
inner regions of the galaxy to be about SxIO^Mq, and the average particle density to be of order 
10^ cm^'^ and the kinetic temperature of molecular gas ~ 50 — 100 K (e.g. Downes & Solomon, 
1998; see Sanders Mirabel, 1996, for a review). The gas in the central ~kpc region of the galaxy is 
likely to have come from distances of the order of 10 kpc when it lost some of its angular momentum 
due to gravitational tidal torques (Barnes & Hernquist, 1992, and references therein; see Barnes, 
2002, for a more recent numerical simulation). 
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The inner kiloparsec region of most young massive galaxies is likely composed of a gaseous disk 
with a mass of several hundred million solar masses, ULIGal being at the extreme end of the mass 
distribution. This gaseous disk is expected to host star formation at a large rate. Some of these 
stars will explode and give rise to shock waves in the gaseous disk which will spawn both more 
star formation and accretion of gas toward the center of the galaxy. We consider these processes 
analytically in some detail in this paper, paying special attention to the effect they might have on 
the evolution of the central parts of the galaxy and on the growth of a central black hole. 

There exists a large body of work on the subject of galaxy mergers, star formation and black 
hole growth e.g. Sanders et al. (1988), Kauffmann k, Haehnelt (2000), Kawakatu & Umemura 
(2002), Granato et al. (2004), Croton et al. (2006), Kauffmann & Heckman (2009), Chen et 
al. (2009), (pi. see Kormcndy & Kennicutt, 2004, for a review), and sophisticated numerical 
simulations e.g. Barnes & Hernquist (1991, 1996), Mihos & Hernquist (1996), Di Matteo et al. 
(2005), Springel et al. (2005), Hopkins et al. (2005), Hopkins & Hernquist (2008). What is 
different in the present work is that we try to capture some of the basic properties of this complex 
system using analytic results for supernova (SNa) remnant evolution and other simple physical 
scalings which are hard to capture in numerical simulations due to the large ratio of galaxy size 
and SNa shell radius. 

The physical system we consider is described in §2 along with the effect SNe have on accretion. 
Bulge formation as a product of SN-induced star formation in the gaseous disk is discussed in §3. 
The main conclusions and uncertainties of this study can be found in §4. 

2. Supernovae-induced accretion in a gas disk 

Numerical simulations of gaseous disks (e.g. Wada & Norman, 2001) find the medium to be 
multi-phase and highly filamentary as a result of star formation and stellar explosion. For the 
analytical calculations in this work, where our primary interest is in average disk properties, we 
consider a simplified disk structure that ignores its filamentary density structure; accordingly, the 
gas distribution is taken to be a smooth function of distance from the center. Many of the results 
reported in this work, as we shall see, have a weak dependence of the interstellar medium (ISM) 
density and therefore the error introduced by the assumption of smooth density field in the disk 
should not be large. 

We consider a disk, roughly 1 kpc in radius, consisting of >10® Mq in gas that came from 
larger radius (^ 10s kpc) due to e.g., tidal interaction with another galaxy. The mean gas density 
in the disk at 1 kpc is ~ 10^ cm~^. Stars form and die in the disk at a certain rate; we take 
the SNa rate in the disk to be dfsn{r)/dA per unit area per year. The interaction between the 
remnant of a supernova with the gaseous disk carves out a cavity in the disk. The gas swept up 
by a supernova is compressed into a thin shell during the snowplow phase, which starts when the 
thermal cooling time scale of the shocked gas is less than the age of the remnant. This means 
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that a certain amount of gas is pushed closer to the galactic center by the supernova shock wave 
and, of course, a certain amount is pushed outward. However, unless gas pushed inwards by the 
supernova shock loses angular momentum it will be pushed back out to a larger radius when the 
shock weakens. We shall estimate later in this section the loss of angular momentum for gas pushed 
closer to the galactic center by the SNa, in order to determine if it is sufHcient to keep the swept up 
gas at a smaller radius. But first, we estimate the amount of gas swept up and pushed to smaller 
radius by a supernova. 

Let us consider a supernova going off in the disk at radius r. The number density of particles in 
the disk at this radius is n(r), the vertical scale height is H{r), and the mean rotation speed of gas 
on its circular orbit is Vgrbir)- During the adiabatic expansion of the SNa shell - the Sedov-Taylor 
phase - the radius, speed and temperature of the shock front are given by 

Rsit) = 3.uf n-^/^^^sf PC (1) 

Vs{t) = 12St^^/^n~^^^El(^ kms-\ (2) 

T{t) = 2.1 X 10^*4^/^71^^/^ ^5/^ K°, (3) 

where t4 = t/lO^yrs, E^i = E/W^^eig, and = n/10^cm~^ is the particle number density 
of gas in the disk at radius r. The Sedov-Taylor phase ends when the radiative cooling time 
tcool = 0.67kT/{nA) is equal to the dynamical time; A ^ 10~^^r~^ergs cm^ s~^ (Blondin et al. 
1998). The radiative phase begins at time 

tsnow ~ 72371-^/^^^5/'^ yrs. (4) 

Subsequently, during the snowplow phase {t > tsnow), ttie evolution is described by (cf. Chevalier, 
1974; eq. 26) 

Rs{t) « 0.8{t/tsnowf^W^^'Eli^^PC^ (5) 
Vsit) 595{t/tsnow)-''-''''nl/^^El(^^kms-\ (6) 
M,{t) « lS3{t/tsnovjfn-'^''Ell/''MQ, (7) 

where Mg is the mass of gas swept up by the SNa remnant. These equations are valid only as long 
as the shell radius is less than the vertical scale-height H: 

^ ,m(cuvr,m ^ .VMCUV?)^^ ^ 3/. / mm ^ -"^ 

^ Vorb V P ; 3 \^iokms-V VlOSMoy ^ ^ 

where rs = r/10^ pc, Cg ~ 1km s^^ is the sound speed, Vt is the RMS turbulence velocity in the 
gaseous disk (produced by SNa explosions and winds from early type stars), and M(r) is the total 
mass enclosed within radius r. 

The Toomre Q parameter for the gravitational instability of the gaseous disk is: 



- 4 - 



where S is the mass density per unit area, Vt,Q = V^/lO^cm s~"^, and Afg = M(r)/10^MQ. We see 
from the above equation that the disk is gravitationally unstable, and wiU support an on-going star 
formation activity. 

Ignoring density gradients in the disk (Rg <C r & i?^ < H), the expansion of a SNa sheh is 
nearly spherically symmetric until the Coriolis or centrifugal force per unit mass becomes of order 
the deceleration of SNa remnant. This occurs when YsQorb ~ d \s/dt, or t ~ 0.69/Qorbj and defines 
a characteristic time when the shell is no longer spherical. 

tshear ~ 7 X lOVa^-fe^s (10) 

where K)rb,2 = Vorh I 100 km s~^. In fact, at this time the magnitude of the shear velocity across the 
shell - Rg dQorb/dlnr - is of order the speed of the SNa remnant. The SNa shell is thus dispersed 
by the shear flow in the disk and mixed with the ambient ISM on this time scale. The radius, 
velocity and mass of the remnant at this time is 

T? f-h \ 1/1 ^-0-25 7710.22 n.3lT/-0.31 

T/ r+ ^ 1 1 »i-0-25 TTn0.22T/0.69^-0.69 „-l /iQ^ 

Vs[tshear) ~ i.irig iig^ Vorb,2^3 km S , [12) 

Msitshear) ~ 7 X 10^ r^-^V-Jf ^i^Of^ (13) 

Note that Rs{tshem-)~H , and therefore the SNa shell is confined within the disk unless r<10pc. 
Since the shell velocity is small compared with the orbital speed when Rg ^ H pc, the swept up 
gas in the disk cannot escape the galactic potential. 

During the snowplow phase the swept-up gas is compressed into a thin shell, forming a hollow 
sphere. Therefore, half of the swept up ISM gas, of mass Ms/2, is pushed closer to the center of 
the galaxy by a distance ~ Rs/2. For an ensemble of SNe going off in the disk the net amount of 
gas that is accreted at r depends on the SNe rate as a function of distance from the galactic center 
(r) and is given by 

Mace - Rs4- f^[2vrri?s]%l ~ M,{RjrfU{r) (14) 



dr \ 2 ' ' dA ^ 

where fsni^) is the cumulative SNa rate within the radius r. We assumed that {MgRsfsn/f) is 
an increasing function of r in deriving the second part of the above equation; otherwise, SNe 
would lead to a net outfiow of gas in the disk to larger distances. Since MgRg cx n^-^^ during 
the snowplow phase (combine eqs. [S] & [3 & substituting for tgnow from eq. ED as long as fsn{i') 
increases with distance faster than there is a net mass accretion even when SNa remnants are 
spherically symmetric. For a Mestel disk the mass enclosed inside radius r increases linearly with 
r, and in that case there is no net accretion - for spherical SNe - if the rate of stellar explosions 
is linearly proportional to gas mass. We shall see below that the condition on fsni^) is relaxed 
when we consider the distortion of SNe remnants by Coriolis and centrifugal forces. Equation (|14p 
for accretion rate is also modified when SNe shells undergo collision before tshear j tbis is discussed 
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below. However, in any case we need to make sure that gas pushed inward by a SNa loses specific 
angular momentum; otherwise, it would be pushed back out when the shock becomes sufficiently 
weak. 

It should be noted that half of a SNa shell has negative angular momentum (as seen by an 
observer at the center of the explosion comoving with the disk), and the other half, with progradc 
velocity field, has positive angular momentum. The magnitude of the total positive/negative angu- 
lar momentum grows during the adiabatic expansion as \L±\ « MsVsr/2 (x l/Vg (the second part 
of this relation follows from energy conservation during the adiabatic expansion phase). However, 
the angular momentum does not increase much during the snowplow phase when MgYg is approx- 
imately constant. The magnitude of the negative angular momentum carried inward by the lower 
half of a SNa shell (|i^ - 1) is quite large and that can lead - as we shall see shortly - to an accretion 
rate in a gas rich disk (like ULIGal) of order a few solar masses per year. 

It can be shown that, for a freely expanding shell, particles on half of the shell with negative 
angular momentum - that have retrograde velocity as seen by a comoving observer at the center of 
explosion - will on average descend a distance of Ar = {SVorb/T^Vg — l)ii^/4r in the radial direction 
(as a result of Coriolis and centrifugal forces in the comoving frame), as long as Vs^Rsld^orb/dlnr]. 
Particles on the other half of the shell with positive angular momentum, or prograde velocity, will 
move outward the same distance. Thus, the shell is continuously distorted with time as a result of 
tidal stretching and Coriolis force. 

Part of the lower half of the shell (lying closer to the galactic center) with transverse velocity 
component in the direction of the orbital velocity has a positive angular momentum with respect 
to the center of explosion, and the part with transverse velocity opposite to the orbital motion has 
negative specific angular momentum. The net amount of negative angular momentum carried by 
the distorted lower half of the shell, during the snowplow phase, as seen by an observer comoving 
with the center of explosion, is 



When the lower half of the shell is mixed - due to shear stretching or collision with other shells 
- the net mean specific angular momentum of the mixed fluid is smaller by ~ 2RsVorb/T^ than a 
particle at the center of the explosion orbiting the galaxy. This is a sufficient amount of negative 
angular momentum to keep the mixed lower half of shell moving on a circular orbit at a smaller 
radius of (r — Rs/2). Thus, the constraint discussed earlier on fsn for inward mass accretion is 
relaxed because there is a net outward transport of positive angular momentum associated with 
each SNa shell. 

The rate of angular momentum transported by an ensemble of SNe in the disk is 



6L « Ar{L_/Rs) « M,i?,K,(,/7r. 



(15) 



L ~ SL{2TrrRs) 



dfsn 

'dl 



2Ms 



TT 



RsVorbiRsMf, 



(16) 
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and the accretion rate resulting from this outward angular momentum transport is 

L 2Msfsn fRsV 



which is similar in magnitude to that given by equation (|14p . i.e. SNe transport mass inward for a 
larger class of functions for fsn{r) than suggested by the discussion following equation ([H]). 

The accretion rate depends on the shell radius (Rs) at the time when SNe shells collide with 
each other or when the shell is dispersed and mixed due to shear flow in the disk. We consider 
both of these cases below. 

Given a supernova rate of fsn per year within radius r of the disk, the mean separation between 
SNe that occurred within time t is 

SNe shells collide when their size is of order the mean separation i.e. Rs{t) ~ dsn- The mean 
collision time, tcoii, is estimated using equation ([5]), and is given by 

f ~ 1 7 V inS ^0.31 ^-0.62^-0.27^1.23 /1q\ 
tcoii ~ i-7 X iU 713 rg yrs. (lyj 

The ratio of the collision time and the time it takes for a SNa shell to be dispersed due to shear 
velocity in the disk is 

^ ^ 40nf-''fL''rf-'X-;l2Elf- (20) 

tcoll 

The accretion rate when tghear/tcou'^^ is given by (using eqs. [TTl [T3l &: [T7|) 



Mace ~ M s{t shear) fs 



Rs{t 



shear- ) 



n 2 



80 fsr^nf-^\f-'^V^lfEk^ Mo yr-^. (21) 



The condition for a remnant not to collide with others before tghear places a limit on the supernova 
rate of 

/.„<3 X 10-3 ny'v'J^.rr'^ E^r' Y^^' - /.n"- (22) 
The second part of eq. (j2ip is valid only when the SNa rate is less than the rate given in eq. ()22p : 
the accretion rate corresponding to this limiting SNa rate is 

Maee ~ 0.3n°-2^£;r -3-°-°Vr,| Me yr-\ (23) 

The accretion rate has a very weak dependence on fsn when the SNa rate is larger than the rate 
given in equation (I22p . i.e. when SNe shells collide before t shear- The reason is that the SNa 
shell radius at the time of shell collision {tcoii oc f^n'^'^) is Rs oc f^n'^^ ^.nd the shell mass is 
Ms (X Rs (X f^n'^'^y therefore, the accretion rate Mace ~ Ms{Rs/rY fsn oc fsn^- For the case of 
t shear /tcoll > the accretion rate is given by: 

Mace « 0.3nf Vr'^sY^s-"-'' M^yr-\ (24) 
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For Mace to be independent of r the density should scale as fsn{i")~^''^T^^'^- However, for a non- 
equilibrium situation in the early phases of galaxy formation and frequent mergers this equilibrium 
density scaling is not applicable. When a large quantity of gas is deposited within a few kpc 
following a merger event or a tidal encounter, the high rate of star formation and SNe at r ~ Ikpc 
would cause a rapid rate of accretion of gas to smaller radii, Mace ~ Ms{tshear)fsn ~ IO^Mq /s„ 
yr~^ after a lag of 10^ yrs, which will continue until star formation and SNa explosions at smaller 
radii start to inhibit this large accretion rate; The subsequent accretion rate would settle down to 
the value given by equation (p^ . 

Note that the cumulative effect of SNe is to create a random velocity field in the disk, but that 
does not automatically ensure accretion. We must have an outward angular momentum transport 
in order for accretion to proceed. An interesting example is that of convective instability in a 
disk; Ryu & Goodman (1992) have shown that disk-convection transports angular momentum 
inward, and therefore the turbulent velocity field associated with it does not lead to any accretion. 
Similarly, the random velocity field in a disk is large when the SNa rate is high and yet because 
shells collide before they are significantly deformed the outward transport of angular momentum 
increases very weakly with fsn (eq. The accretion rate in terms of an effective Q-viscosity, 

u = aRs(tcoii)Vs{tcoii), in the limit that fsn > fsn'' is given by 

Mace ~ 2TTm.pnHu ~ 1.2a(i//0.1r)n0-39/r'^5i'^^3 Moyr"!. (25) 

Comparing this with equation (f24l) we see that a ~ 1 for fsn ~ fsn' (given by eq. [22]) . and a is 
smaller for a larger SNa rate, although the effective a is dependent on r. 

If SNe remnants punch through the disk in the vertical direction, but are still confined by 
the galactic potential, then some fraction of gas leaving the disk will eventually fall back onto the 
disk at a smaller radius and contribute to the net accretion rate. A SNa is confined to the galaxy 
provided that 

Rsimm{tcoll,tshear}) < H OT Vs{Rs = H) < Vorb- (26) 

If many SNe shells collide and coalesce before they are dispersed they would form super-shells 
and if their velocity is sufficiently high they can escape the galactic potential. Otherwise, these 
super-shells would also contribute to transporting gas to smaller radius. 

3. Star and bulge formation and growth of a central black hole 

The possibility of SNa-induced star formation has been discussed in numerous contexts, from 
the early universe to the present day Milky Way (e.g. Woodward 1976; Bedogni & Woodward 1990; 
Yamada & Nishi 1998; Mackey et al. 2003; Bratsolis et al. 2004; Joung & Mac Low 2006; Johnson 
& Bromm 2006; Sakuma & Susa 2009; Leao et al. 2009; Nagakura et al. 2009). Here we consider 
how this process may compete with the fueling of black holes and contribute to the formation of 
galactic bulges. 
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For a Miller-Scalo initial-mass function (IMF) for stars (Scalo, 1986): 

1.9(M/O.OlM0)-"i O.OIM© < M < O.OSMq 

(M/O.O8M0)-°2 O.O8M0 < M < 0.5Mq (27) 

6.25-"2(7if/o.5M0)-"3 M>0.5Mq, 

with ai = 0.3, 02 = 1-8 ± 0.5, = 2.3 it 0.7 (the parameters are taken from Kroupa, 2001), and 
A^* = / dM dN/dM. The mass fraction in high mass stars {M^'>^Mq), capable of SNa explosion, to 
the total star mass is about 0.3, and the fraction by number is about 4 x 10^^. Thus the expected 
SNa rate is 0.3 yr^-*^ for the Miller-Scalo IMF and a star formation rate (SFR) of 10 M0 yr^-*^, 
which is of the order observed in ULIGals and needed for forming galactic bulges in L* galaxies. 

Indeed, this is consistent with the following simple estimate of the star formation rate in the 
gaseous disk: 



dN 
dM 



4.5iV* X < 



SFR ~ /^Mdisk/tff, (28) 

where Afdisk is the disk gas mass, tg is the free-fall time, and /* is the efficiency with which gas is 
turned into stars on a free-fall time, taken to be of the order of /* ~ 0.01 (e.g. Krumholz & Tan 
2007). Then, for a disk of 10^ - 10^ Mq, the SFR is about 1-10 Mq per year, assuming a density 
of 10^ cm^ 



3 



If a similar fraction (~ 1%) of gas in SNe shells is turned into stars (before shells collide) then 
the resulting star formation rate would be ~ lAO Mqu'^'^'^ E^fr^'^ fg^ yr"^, and that would result 
in a SNa rate of ~ 5n!^'^'^ E^fr^'^ f^^ yr"""^ if the IMF for these "daughter" stars were given by 
equation (f27|) . However, we show below that formation of massive stars - those capable of SNa 
explosion - is suppressed in SNe remnants when /s„>0.5. 

Star formation in a SNa shell is suppressed on large length scales due to the transverse relative 
velocity gradient in the remnant. We calculate this length scale as well as the Jean's length, and 
estimate the rate of star formation in SNe-remnants. 

The velocity field in a SNa remnant seen by a comoving observer in her neighborhood is: 
6Vs ~ {Vs/Rs)Sr; where 6r is the position vector, tangential to the shock front, pointing from the 
comoving observer to a point in her neighborhood. For a gas clump to be able to collapse, the 
differential velocity across the clump, \SVs\, should be smaller than the gravitational escape speed, 
i.e. for a clump of size £, \6Vs\ ~ Vs{i/Rs) < [GMs{l/ Rsf /UY^'^ or 



Using equations ©-([T]) we find 



~ 9 X IQ-^^ 



4.9 



_6_ 31 
'17 17 17 



Eil M0, (30) 
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We see from equation ([6]) that at t/tgnow ~ 137n^'^^ E^f^ the shell speed is Vg = 20 km s ^, and at 
that time Mi ~ 0.2n^''^^E^i Mq; the maximum star mass scales as V~'^'^. 

Star formation is disrupted when shells collid^. The average speed of SNe shells when they 
collide, at time tcoii, is Vs ~ 14 km s~^fg^r^^'^n^^'^E^^ (obtained from equations 6 fc [T9]l . For 
a small SNa rate, fsn < 0.5 yr~^, Vg can drop down to a value where massive stars capable of 
SNa explosion can form before SNe shells collide and star formation is disrupted; however, since 
lower mass stars form first in SNe remnants they can perhaps significantly suppress the formation 
of more massive stars. At a higher SNa rate, formation of stars more massive that a few solar 
mass is suppressed; massive stars could still form in the disk after the turbulence generated by shell 
collisions has subsided and gas has cooled down. This suggests that SNe explosions in a gaseous 
disk might occur in waves of high and low activity, and during periods of high activity Miller-Scalo 
IMF is truncated above a few solar mass, and even during periods of lower SNa rate the IMF could 
be more bottom- heavy than the standard IMF. 



3.1. The thermal state of the pre-shock gas and the shock-induced stellar IMF 

A firm lower limit to the shock front speed, which in turn sets an upper limit to the masses of 
the stars which may form behind the shock front, is set by the sound speed of the medium upstream 
of the shock front when SNa shock weakens and turns into a sound wave. In order to calculate the 
sound speed of the upstream gas, we must determine the thermal state of this gas, and we discuss 
two possible cases in this subsection. 



3.1.1. High supernova rate 

For the first case, we consider a galactic disk with a high star formation rate (SFR), in which 
the dominant process affecting the thermal state of the gas is photo-heating of the gas by massive 
stars which eventually explode as SNe. Considering that these stars will each live for roughly 10^ 
yr, the average distance between these massive stars within radius r of the disk, following equation 
(18), is 

dsn~0.6-^/,„(r)-i/2pc. (31) 
Ikpc 

The number of ionizing photons emitted per second by such massive stars is > 10^^ (e.g. Osterbrock 
& Ferland 2006). Taking the average gas density to be ~ 10^ cm^^, we find the radius of the 
Stromgren spheres surrounding massive stars (Stromgren 1939) to be of order 1 pc. Thus, we 
expect that for an average SNa rates of /sn~ 0.4 yr~^ the H II regions surrounding the massive stars 



^Collision between SNe shells disrupts star formation because of turbulence generated in these collisions. Once 
turbulence dies out - in about a shock crossing time - star formation can resume provided that the merged shell is 
not hit again by a high speed shock. 



-lo- 



in the disk will overlap, and therefore SNe shocks will propagate into such photoionized regions. 
Largely independent of the metallicity of the gas, the temperature in such a photoionized region 
will be of the order of 10^ K, and the sound speed ~ 10 km s~^. Therefore, SNe shocks in this case 
will generally dissipate and turn into pressure waves once they have slowed to speeds of ~ 10 km 
s~^. This suggests, following equation (29), that the IMF of stars formed in the material swept up 
in SNa shocks is likely cut-off at a few solar masses. 

3.1.2. Low supernova rate 

We next consider the case of a lower SFR, corresponding to 0.4 yr~^, for which the 

distance between massive stars is greater than their average Stromgren radii. In this case, SNa 
shocks will generally propagate into so-called relic H II regions in which gas that was previously 
photoionized by the progenitor star is recombining and cooling; equations ([5]) & (l6|) show that, 
in general, i?s>lpc (Stromgren radius) when the SNa shell speed has dropped to ~ 10km s~^. 
However, even for larger Stromgren radii, the upstream relic H II region gas can in some cases cool 
to temperatures well below lO^K, which allows for the formation of stars with masses greater than 
a few solar masses. 

The results of a numerical calculation of relic H II region gas temperature ahead of the shock 
are shown in Fig. 1. This calculation assumes that the density of the gas remains constant at the 
fiducial value of n = 10^ cm~^, and that the gas cools only radiatively through atomic transitions of 
metals, which is a reasonable assumption for the case we consider here (but see Jappsen et al. 2009, 
for the case of star formation in a more isolated environment). In principle, the radiation emitted 
by the shocked gas can send a radiative precursor ahead of the shock and heat the upstream gas, 
however, we neglect this effect as at the shock velocities at which the shocks stall, i.e. < 20 km 
s~^, radiative precursors do little to ionize or heat the upstream gas (Shull & McKee, 1979). 

Figure 1 shows results for four different metallicities: 10^^, 10^^, 10^^, and 10~^ Zq. The 
metallicity-dependent cooling rate of the gas is taken from Mashchenko et al. (2008), who provide 
a fitting formula to the cooling function calculated by Bromm et al. (2001). The left panel of Fig. 
1 shows the temperature of the upstream gas as a function of the time from the death of the central 
star. The right panel shows the square of the Mach number for the shock front. The SNa shock 
will stall once this ratio approaches unity - at this point the shock will turn into a pressure wave. 

For metallicities $10~^ Zq the SNa shock stalls at a time igtaii ~ 10^ yr, when the shock 
velocity is still ~ 10 km s~^. However, for higher metallicities the shock stalls at later times when 
Vs has dropped to a smaller value and thus the maximum mass of SNa-induced stars can become 
M/-^ 5 Mq for shocks that last for > 2 x 10^ yr (see eq. 30). We note also that for the case 
we consider here of fsn ^ 0.4, tcoU can easily exceed tstaih and therefore shell collisions will not 
interfere with SNa-induced formation of these more massive stars at late times. 
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Fig. 1. — The temperature of the medium upstream from the shock front (left panel) and the ratio 
of the squares of the upstream sound speed and shock front speed during the snowplow phase (right 
panel), in a cooling relic H II region with gas density of 10^ cm~'^. The shock wave turns into a 
compression wave when the ratio drops to ~ 1, which occurs at t > 10^ yr for each of the gas 
metallicities considered here: the curves from left to right correspond to metallicities of 10^^, 10^^, 
10"'^, and 10~^ Zq, respectively. Note, however, that for lower upstream gas densities this occurs 
at earlier times, due to the lower cooling rate of the relic H II region gas. 



3.2. Stars, bulge and black hole 

As shown in the last Section, for /sn~ 0.4 yr~^, stars forming in SNe remnants have peculiar 
velocities of order 10-20 km s~^ in radial direction with respect to the center of the explosion. 
When viewed from the galactic center the velocities of these newly formed stars in SNa remnants 
would have a random velocity dispersion of ~ 10-20 km s~^, and therefore these stars would tend 
to form a bulge at the center. The velocity dispersion of the newly formed stars would increase with 
time as these stars are subjected to the stochastic gravitational field of the SNa remnant network. 
Moreover, these stars will suffer some hydro dynamical drag on their way out of the gaseous disk, and 
also will be subject to gravitational drag that will modify their velocity dispersion (eg. Artymowicz, 
1994; Nayakshin & Cuadra, 2005); the former is likely a small effect due to the small cross-section 
for star-gas interaction, but the latter can be a significant effect and needs to be included in the 
calculation to determine the true random velocity distribution of stars formed out of SNe remnants. 

The supernova led accretion would also deposit gas in the central parsec region of the galaxy at 
the rate given by equation (j24p : at distances smaller than ~ Ipc the magneto-rotational viscosity 
(Balbus and Hawley, 1991) is expected to be effective in transporting gas to the black hole at 
the center. We note that MRI might also operate at larger radius due to heating and ionization 
produced by SNe. 
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Let us assume that a fraction /* of supernova remnant mass is converted to stars before shells 
collide with another shell. Using equations ([7]) and (|19p we estimate the star formation rate in SNa 
remnants to be 

M, « (1.5 X lO^Mo) /, /,Vn0-54£;0^Vi-i yr"! (32) 

Therefore, the ratio of star formation & accretion rates for /^n^JO.Syr"^ is ^BxW^ fiffg^iT's'^ E^i'^^s'^ 
(eqs. [24l&: [32]) . For /* ~ 10~^ and fsn ~ 1 yr~^ this ratio is of order a few hundred. In the case of 
low SNa rate (when the remnant survives until its velocity drops to ~ 10 km s^^ and then turns 
into a compression wave) the ratio is 3 x lO^/^ng ^£'^^'^ ''r|, which is also of order a few hundred 
for ~ 10~^; this last expression was obtained by using eqs. [5] & [7] for Rg & Ms corresponding to 
shell velocity of 10 km s~^ (eq. ^ to calculate star formation and accretion rate (eq. [T7|) . 

This ratio of star formation rate to accretion rate is similar to the reported ratio of bulge and 
BH mass in galaxies (Gebhardt et al. 2000; Ferrarese &: Merritt 2000), and the scenario we have 
described offers a plausible physical explanation for this correlation. We note that a number of 
well known feed back processes have been left out in the calculations presented in this paper (cf. 
Cattaneo et al. 2009), and these might significantly modify star formation and accretion rates. 

4. Summary 

We have analyzed the effect of supernovae occurring within the central kpc region of a gaseous 
disk on the formation of stars and transport of gas from ~ Ikpc to a few pc of the galactic center. 
The outward transport of angular momentum facilitated by SNe explosions allows for the inward 
transport of gas that feeds the central black hole. This is a process which may take place quite 
generically in any galactic disk hosting star formation, although it may not be the dominant process 
affecting black hole accretion. 

We have shown that associated with the inward transport of gas swept up by SNe is the shock- 
induced formation of stars which are born with a random peculiar velocity of ~ 10 km s~^. This 
velocity dispersion increases with time as a result of the stochastic gravitational field associated with 
filamentary SNa remnants; The stars formed in SNa remnants contribute to the stellar population 
of a central bulge. 

Due to the divergent velocity field of an expanding SNa shell there is a maximum length scale 
for fragmentation of shells or an upper limit to the mass of stars formed in SNa remnants; the 
SNa-induced stellar IMF is cut-off above a few solar masses and more massive stars can only form 
if and when a SNa shell slows down to ~ 10 km s~^. 

We note that numerous observations have suggested connections, such as we have considered in 
the present work, between star formation, black hole accretion, and the formation of a stellar bulge. 
Heckman et al. (2004) suggest that star formation and black hole accretion rates are correlated, 
and Chen et al. (2009) report an empirical relation between supernova rate and gas accretion 
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rate on the central black hole (see also Xu & Wu 2007). Furthermore, Page et al. (2001) report 
observations suggesting that central black holes and stellar spheroids form concurrently, and Genzel 
et al. (2006) describe observations of a galaxy hosting both an accreting black hole and a central 
stellar bulge, with no evidence of a major merger. Hydrodynamic simulations have demonstrated 
that SNa feedback may produce spherical distributions of stars in dwarf galaxies (Stinson et al. 
2009) and in the inner portions of the Galactic bulge (Nakasato & Nomoto 2003) . We would like to 
point out the recent work of Wang et al. (2009) that models SNe induced turbulence as an effective 
viscosity and describes the evolution of a gaseous disk. 

A limitation of this work is that we have ignored radiative feedback effects which are known 
to control the steady state accretion rate onto the black hole (e.g. Ostriker et al. 1976, Proga 
et al 2008, Milosavljevic et al. 2008) and probably also affect the formation rate of stars in the 
central kpc. Ultimately, large scale simulations resolving the long term evolution of individual 
SNe remnants, star formation, the multiphase ISM, and the feedback effects of accretion onto a 
central black hole will be required to more fully elucidate what role SNe-induced accretion and star 
formation play in galaxy formation. 
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